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Many and very general arguments indicate that the event horizon behaves as a stretched mem- 
brane. We explore this analogy by associating the Gregory-Laflamme instability of black strings 
with a classical membrane instability known as Rayleigh-Plateau instability. We show that the key 
features of the black string instability can be reproduced using this viewpoint. In particular, we get 
good agreement for the threshold mode in all dimensions and exact agreement for large spacetime 
dimensionality. The instability timescale is also well described within this model, as well as the di- 
mensionality dependence. It also predicts that general non-axisymmetric perturbations are stable. 
We further argue that the instability of ultra-spinning black holes follows from this model. 
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The existence of black holes is perhaps the most dra- 
matic prediction of Einstein's theory, the very concept of 
which makes full use of the non-linearity of the equations 
and also of our notion of space and time. Despite being 
apparently very complex objects, black holes can be as- 
sociated with many of the familiar quantities of everyday 
physics. The first major breakthrough in this direction 
was hinted at by Bekenstein [lj, who conjectured that 
black holes are endowed with thermodynamic properties, 
namely with an entropy proportional to its area. That 
black holes are thermodynamics entities was established 
once and for all by Hawking Q, who verified explicitly 
that black holes radiate, and therefore have an associ- 
ated temperature. This analogy carries over to higher 
dimensional scenarios, which seem to be a prerequisite 
for consistency in many modern theories. In the general 
higher dimensional gravity theories, there are other ob- 
jects called black branes with event horizons. These 
are basically extended black holes: the horizon, instead 
of having the topology of a sphere, can have for instance 
the topology of sphere times a line - a cylinder. The anal- 
ogy with thermodynamics can still be formulated, and in 
fact one can even go a further and formulate an anal- 
ogy with hydrodynamics 0]. The thermodynamic de- 
scription, both for black holes and black branes, is based 
on the four laws of black hole dynamics formulated by 
Bardeen, Carter and Hawking || . The first law (we will 
take for simplicity uncharged, static objects) describes 
how a black hole, characterized by its mass M, horizon 
area A, and T = l/(327rM), evolves when we throw an 



infinitesimal amount of matter into it: 
dM = TdA . 



(1) 
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The second law states that in any classical process the 
horizon area must increase, dA > 0. It is very tempting 
to associate these two laws with the first and second laws 
of thermodynamics, respectively, in which case T would 
be proportional to a temperature and A to an entropy 
(the other two laws of black hole mechanics also have a 
correspondence with the zero and fourth thermodynamic 
laws). The final ingredient to proceed consistently with 
this association was given by Hawking |2|, who realized 
that black holes are indeed radiating objects and that one 
can indeed associate them with a temperature Th = 4T. 

We can also argue that Equation 0} can be looked at 
as a law for fluids, with T being an effective surface ten- 
sion |6|]. Regarding the event horizon as a kind of fluid 
membrane is a position adopted in the past Q. The 
first works in black hole mechanics actually considered 
T as a surface tension (see the work by Smarr [j| and 
references therein), which is rather intuitive: in fluids 
the potential energy, associated with the storage of en- 
ergy at the surface, is indeed proportional to the area. 
Later, Thorne and co-workers [jj developed the "mem- 
brane paradigm" for black holes, envisioning the event 
horizon as a kind of membrane with well-defined me- 
chanical, electrical and magnetic properties. Not only 
is this a simple picture of a black hole, it is also use- 
ful for calculations and understanding what black holes 
are really like. There are other instances where a mem- 
brane behavior seems evident: Eardley and Giddings |IJ, 
studying high-energy black hole collisions found a soap 
bubble- like law for the process, while many modern inter- 
pretations for black hole entropy and gravity "freeze" the 
degrees of freedom in a lower dimensional space, in what 
is known as holography In .4] it was shown that a 
"membrane" approach works surprisingly well, yielding 
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precisely the same results as the AdS/CFT correspon- 
dence. There have also been attempts to work the other 
way round: computing liquid surface tension from the 
(analog) black hole entropy 

We will consider here an interesting class of black ob- 
jects without spherical topology, the black branes, and 
treat them as fluids with a surface tension subjected to 
JU and without gravity. A broad class of these objects 
are unstable against gravity, in a mechanism known as 
Gregory-Laflamme instability ^3] • This is a gravitational 
instability along the "extended" dimension, pictured in 
Figure n We will show that we can mimic most aspects 
of this instability using the fluid analogy. 

Take a black string with radius Ro, and extended along 
a direction z. Then the Gregory-Laflamme mechanism 
makes any small perturbation with wavelength A of the 
order of, or larger than, the radius of the cylinder Rq 
grow exponentially with time. This is a very robust long- 
wavelength instability. The details for the instability can 
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FIG. 1: Black strings and fluid cylinders are unstable to per- 
turbations on the extended dimension, i.e., along the axis of 
the cylinder. Ripples propagating along this axis grow expo- 
nentially with time for wavelengths of order of the radius of 
the cylinder. 

be found in the original papers [f2T |. For wavelengths 
larger than a threshold A c , the instability appears. In the 
general setup of black strings with D— 1 spatial directions 
with radius Rq and a transverse direction z, one finds 0] 
that the threshold wavenumber increases with dimension 
number D and so does the maximum growth timescale. 
For very large number of dimensions, the threshold mode 
behaves as Il3l 



k c = 2ir/\ c ~ VD/R 



(2) 



Quite remarkably, we can find a similar instability for 
fluids with surface tension, as shown in Plateau's cele- 
brated study 0| on the stability of bodies under the 
influence of surface tension. He established a fundamen- 
tal result of classical continuum mechanics: a cylinder 
longer than its circumference is energetically unstable 
to breakup. This result was put on a firmer basis by 
Lord Rayleigh ^{|, who computed the exact instability 
timescale for the problem. The reasoning for the appear- 
ance of an instability is the following: consider a small 
disturbance of a long cylinder (we take its axis as the 
z— axis) of fluid with radius Rq and height z. Consider- 



ing a small axisymmetric perturbation along the surface 
of the cylinder, we write for the disturbed cylinder 



r(z) — R Q + e Ri cos(fcz) + e 2 R 2 , 



(3) 



where e measures the perturbed quantities (R 2 is a sec- 
ond order quantity, and its usefulness will be understood 
shortly). The volume of this cylinder can be easily com- 
puted to be 
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2RoR 2 
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(4) 



R 2 

If we impose constant density one must have R 2 — — . 
With this condition, the surface area of the disturbed 
cylinder is 



A = ZTT 



2Rn 



e 2 R 2 
2Rq 



(k 2 R 2 1) 



The potential energy per unit length is therefore 
ire 2 R 2 



P = 



2R 



{k 2 Rl - 1) T. 



(5) 



(6) 



We conclude that the system is unstable for k < 1/Rq, 
since in this case a perturbation of the form Q decreases 
the potential energy. 

One can further show that non-axisymmetric perturba- 
tions, with profile r(z, <fi) = Ro + eRi cos(kz) cos(m4>) + 
e 2 i?2 (where m is an integer that identifies the angular 
mode) are stable for any m ^ ,15] . The reason is that 
the potential energy for these non-axisymmetric modes 

is given by P — ^^r!, 1 (k 2 Ro — 1 + ti 2 ) T, which never 
decreases for m ^ 0. 

We can generalize the Rayleigh-Plateau construction 
to a general number of dimensions. Take a hyper-cylinder 
with D — 1 spatial directions with radius Ro and a trans- 
verse direction z (for the previous example D = 3). The 
axisymmetric threshold wavenumber is 



R k c = VD - 2 with k c = 2ir/\ c . 



(J) 



For wavenumbers smaller (larger wavelengths) than this 
critical value the cylinder is unstable. Moreover, as in 
the original Rayleigh-Plateau situation, only symmetric 
modes seem to be unstable; non-axisymmetric modes are 
in general stable. Therefore, the Rayleigh-Plateau in- 
stability, like the Gregory-Laflamme instability, should 
disappear for modes other than the s— modes. This 
was recently conjectured by Hovdebo and Myers 
using an argument based on the relation between the 
thermodynamic and the Gregory-Laflamme instabilities 
|l7| . Kudoh [l8| has explicitly verified that the Gregory- 
Laflamme instability only affects s— modes. 

To motivate quantitatively the suggested associa- 
tion between the Rayleigh-Plateau and the Gregory- 
Laflamme instabilities, it is important to compare the 
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dependence of the threshold wavenumber kRo on dimen- 
sion D, for both instabilities. This is done by comparing 
© with 0). In Tabic U we list the value of the Rayleigh- 
Plateau threshold mode Rok c for several dimensions D. 
We also list the threshold wavenumber for the Gregory- 
Laflamme instability, with values taken from 0] • There 
is good agreement between them. In the large Z3-limit 
there is exact agreement: both the Gregory-Laflammc 
and the Raylcigh-Plateau critical wavenumber behave as 
k c Ro ~ \f~D. This is, we think, a non-trivial check on the 
conjecture that black branes behave as fluid membranes 
with surface tension. We can further compare the evo- 
lution of the instability timescale with its wavelength, 
and study the dependence of the instability timescale 
on the spacetime dimension. We compute the Rayleigh- 
Plateau instability timescale using fluid dynamics, fol- 
lowing Rayleigh [15j. For a 3-dimcnsional cylinder, and 
assuming the perturbation goes as R\ ~ e ot , Rayleigh 
gets the following expression for fi: 



9. 2 



T ikRoJ' {ikRo) 
pRl J (ikR ) 



[l - k 2 R^ 



(8) 



with J being a Bessel function. We have generalized 
Rayleigh's procedure for higher dimensions, and the re- 
sults are shown in FigureEl We assumed an effective sur- 
face tension and density associated to the temperature 
and energy density of a Z?-dimensional Schwarzschild 
black hole. The qualitative behavior of the instability 
timescale matches surprisingly well that of the Gregory- 
Laflammc 12] one. In particular note that (i) the max- 
imum growth rate grows with the number of spacetime 
dimensions, and (ii) the corresponding wavenumber also 
grows with dimension number D. Moreover, the lo- 
cation of the threshold wavenumber is exactly as pre- 
dicted with the energy argument (see eq. ©)■ There is 
one discrepancy only: the maximum instability for the 
Rayleigh-Plateau case, i.e., the maximum Q is approxi- 
mately one order of magnitude larger than the maximum 
of the Gregory-Lafiamme. This could be due to the com- 
plete neglect of gravity in the outside of the cylinder (in 
the Rayleigh-Plateau analogy). Indeed if one included 
gravity effects, a redshift was bound to occur, thereby 
lowering f2 |igj . 

It is interesting to ask what is the endpoint of the 
Rayleigh-Plateau instability. In four-dimensions, the 
Rayleigh-Plateau instability makes a cylinder pinch off, 
giving rise to drop formation (see for instance the original 
works 0, 0] or the review work by Eggers [2£| ) . To see 
that this is expected, consider af = 3 cylinder breaking 
up at a Rayleigh-Plateau length, i.e., a cylinder of height 
z = A c = 2ttRq . Can this give rise to a drop of the same 
volume? Same volume implies that the radius of the final 
spherical drop would be 1.6765i?o- This makes the sur- 
face area of the drop smaller than the corresponding for 
the cylinder, Sdxop — 0.8947S , cy i, and thus it makes sense 
that a drop is the final state. Now, if we play this geomet- 




FIG. 2: Rayleigh-Plateau instability of a hyper-cylinder in 
several dimensions. Here the effective surface tension and 
density were chosen to match those of a higher-dimensional 
non-rotating black hole. 



rical game in higher dimensions, we conclude that there 
is critical dimension between D = 10 and D = 11. For 
D > 11 a sphere seems not to be the favored endpoint, 
since it no longer has less surface area. Once again, there 
is a very similar phenomena in the gravity case (see ,e.g., 
the discussion in [lil EH Ei| ) . Take a black string break- 
ing up at a Gregory-Lafiamme length. Sorkin found a 
critical spacetime dimension (d = D + 1) between d = 13 
and d = 14, above which the black string is no longer en- 
tropically unstable against the formation of a spherical 
black hole 21]. It would be interesting to further study 
this issue. In particular, it is important to understand 
what is the endpoint of the Rayleigh-Plateau instability 
for D > 11, and to find if in the fluid model there is a 
new branch of solutions that would be the analogue of 
the non- uniform black string solutions of 23], and if so 
to study their stability. 

It has been shown by Horowitz and Maeda [3] that 
pinch-off, if it occurs at all in the black string case, must 
do so in infinite affine time. This immediately suggests 
that an attractive endpoint could be non-uniform black 
strings H3,0|- Now, the breakup of liquid jets with sur- 
face tension in the absence of viscosity is known to hap- 
pen in finite time, but the inclusion of viscosity (which 
seems a necessary ingredient to model realistic black ob- 
jects [2^ ) may change this [2(|, so even this unexpected 
feature might be discussed within this analogue model. 
It is quite amusing that some of the dynamical features 
of the instability of black strings have already been ob- 
served in liquids: the final state of some liquid bridges 
(finite-size liquid cylinders), unstable under the Raylcigh- 
Plateau instability, is a non-symmetric state (see |2(j and 
its references) and the breakup of liquid jets is quite gen- 
erally a self-similar phenomena 20]. 

The interpretation of T as surface tension can also im- 
prove - and strengthen - our understanding of the in- 
stability of ultra-spinning black holes (Myers-Perry |27j 
black holes with high rotation), a conjecture recently 
made by Emparan and Myers in connection with the 
Gregory-Lafiamme instability [2^]. Take a slowly rotat- 
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TABLE I: Dimensionless threshold wavenumber kRo for the Rayleigh-Plateau instability of a higher dimensional fluid cylinder, 
and the corresponding threshold wavenumber for the Gregory-Laflamme instability (data taken from [X^])- 



D spatial dimensions 


4 


5 


6 


7 


8 


9 


49 


99 


Rayleigh-Plateau 


1.41 


1.73 


2.00 


2.24 


2.45 


2.66 


6.78 


9.80 


Gregory-Laflamme 


0.876 


1.27 


1.58 


1.85 


2.09 


2.30 


6.72 


9.75 



ing black hole. As the rotation rate increases, the surface 
becomes flatter at the poles until zero Gaussian curvature 
occurs. But a liquid drop (at least in four dimensions) 
develops instabilities before zero Gaussian curvature is 
reached [13 ■ Assuming a correspondence between rotat- 
ing liquid drops and rotating black holes, one may well 
expect ultra-spinning black holes to be unstable. We 
note that this very same reasoning was applied by Smarr 
many years ago. In four dimensions, it seems that 
the upper Kerr bound in the angular momentum is small 
enough to avoid the development of such instabilities. 
However, in dimensions higher than six, rotating black 
holes have no Kerr-like bound, and the instability might 
well set in. 

It is possible that an event horizon behaves dynam- 
ically much as a fluid interface without gravity, as we 
have shown. If this is so, the second law of black hole 
dynamics should be something like a soap-bubble law, 
sphericity is preferred for it is the minimum energy shape 
(which would justify why most stable solutions in GR 
have spherical topology) . Finally, there are many ways to 
extend these results, by including additional effects, like 
charge or rotation in the problem. One would naively ex- 
pect either of these to increase the instability, but a more 
careful study has to be done. Take rotation for instance. 
The centrifugal force, scaling with r, contributes with a 
de-stabilizing effect, since there is an increased pressure 
under a crest but a reduced pressure under a trough. So 
in principle, the threshold wavelength should be smaller 
However, there is no simple reasoning (as far as we 
know) to tell the effects of rotation on the strength of 
the perturbation, i.e., on the instability timescale. This 
is partly because the effective surface tension of black 
hole itself depends on rotation, decreasing with increas- 
ing rotation. This would be an interesting direction for 
further study. 
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